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I. Il^TBODUCTION 

1 

. . I 

Centaurus X-3 was the first discovered X-ray source whose intensity is 
pulsed and eclipsed (Giacconi ^ 1971 1 Schreler 1972). Though the 

X-ray properties of the Cen X-3 system have been intensively studied, the optical 

properties remained unstudied until Krzeminski’e (1973) identification of the 

I 

optical counterpart of Cen X-3. Observations of the X-ray eclipse duration, 
pulsation fre<3,uency variation, and pulse ELrrival time variation have been used 

pature of the system and, to predict the nature of the 
optical companion e^id. X-ray source (Schreier ^ al . 1972, Sofia 1972^ Wilson 
1972, Os^i4972, .pavidson and Ostriker 1972, and M 9 C^uskey and Kondo 1972). 

It is purpose, qf^, 'll, hi , paper to present an observational determination of 
several system parameters based on a prellmlneury interpretation of tjhe light 

i ' 

curve of j Xr zeminsk^ * a star . 

I 

/ . i 

\ ' Hi OBSERVATIONS REDUCTIONS 

■ I 

I 

The observations presented here were obtained with the Cerro Tololo 

, I 

Inter-American Observatory 1 meter (Yale) telescope and the 6l cm. (Lowell) 
telescope In the period 25 January 197^ ito 12 February 197i*. A 1P21 (CTIO 

V ■ * • 1 =■■.•;'■■■ . • 

no. 57) was used in pulse counting mode on both telescopes with the V filter 
of the CTIO no. 2 UBV filter set. The standard CTIO photometers were used, 
photometer no. 2 being used on the 1 meter telescope and photometer no. 1 
on the 6l cm. telescope. ^ 

Observations at the telescope were made in the order: sky-comparison- 

Krzeminaki's star, being sure to bracket all observations of Krzeminski's star 
with comparison observations bracket all comparison measures by sky 

measures. All obs^i^yations were one minute integrations. Observations were 
made in this se<iuencc until approximately 10 measures each of the sky, Krzeminskl * s 

I ' • ■ 

; / 
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st€kT, and the comparison were obtained.' The primary comparison observed was 
no. 9 of Brucato, iCristlan» and Westphal (1972) and the secondary standard was 
no. 6 of the same paper. 

t 

In the data reduction sky subtraction was effected by taking the mean of 
adjacent sky measures and applying this to the bracketed measures of the com- ^ 
parison and of Krzeminski’s star. The differential magnitude of Krzeminski’s 
st^ was then fovipd. wltjji respoct to the mean of the^ comparison star measures 
bracketing; th^, measure, of Krzeminskl *8 star. . 

Two problems were encountered in these procedures. The first was due to 

the crowded field, in.. w|iicj> IQrzeminski's star la located. Observations obtained 

with a 20 arc second aperture differed systematically, from those obtained with 

a 15 arc secopd apertiire. A special sequence of observations was made on two 

nights with the 1 me;j^^r telescope fiuid, combined with a previous night’s observations, 

* 

when both apertures on the 1 me'jjer^ telescope, had been employed to determine 
the aperture effect bp 0.056 ± 0.005 mag. This correction has been applied 
to all observations obtained with the larger apertures. The size of this 
effect corresponds to the presence of a star of nagnitude l6.5 in the larger 
apertures. Residuals from the light curves described in the following section 
show no statistically significant difference between the two sets of observa- 
tions, af;ter co?:,3;pction of the larger aperture set. , 

The ppcondj prpbjLpm, was that the sky subtraction method failed for obser- 
vations made h|epr, moonrise and moonset . . , !Tp overcome this, we interpolated graphi- 
cally in the Ipgnriifc.lm of the sky brightness to determine the value of the sky 
background at the time, of each observation. The difference between the logarithmic 
graphical procedure and our standard jbec^ique was as large as 0.03 magnitude. 




Haying taken these effects into account the mean and Internet, error of 

j ■ j 

each of the 40 minute runs was computed* The standard deviation of an ohser- 
vation of unit weight is 0*008 magnitude. Figure 1 presents these means from 

which it may he seen that the amplitude of the light curve is approximately 

! 

0.08 magnitude, in agreement with the initial announcement hy Krzeminski (1973) 
hut in disagreement with the value of 0*12 magnitude reported recently hy Mauder 
( 1974 ) and Krzemsixikl (197^*)*' However, if the extreme veJ.ues of the present 

j 

observations are employed in determining the amplitude, then a value approaching 
0.12 magnitude is found. . ..j, 

j 

III. LI6IIT' CURVE SYHTHESirS AND SYSTEM PARAMETERS 
We have determined the mass ratio' q = and the inclination. 



for synthetic light curves over a range 6f (q, sin i) values. These light curves 

■i . ■ - . . [ 

are computed from formulas given hy Kopal and Kitamura ( 1968 ) for the light 
changes of close hineiry systems. The model employed is the standard Roche 

^ i i i V i 

model wherein the steurs are assumed to I he of infinite central condensation, 

' i 

rotating synchronously, and to have the rotational and orbital angular momenta 
parallel, ^e observed surface brightness is taken to obey von Zeipel’s theorem 
and a n^inear limb darkening law. We have assumed that the primary star fills 

i 

its Roche lob^, ,„.Th! 5 . effect of this assumption is that the mass ratios and inclin- 
ations determined will be minimum values. 

We have eijiployed the. NLTE, line blanketed, plane parallel,, hydrostatic 
model atmospheres of Mihalas (1972) to' compute the gravity darkening coefficients 
Tq and values of which, depend upon the spectral type of the primary. 



The following spectral types for Krzemlnski ' s star have been suggested: 

09 - 09.5 V, BO I-III, and BO II by Krzemlnski (197^) baaed on the UBV colors 
of the star; 09 III-V (Vidal et 197^); 09-5-B0.5 Ib by Rickard (197^) 
from, spectra; and 06.5 III-V (Osmer, Hiltner, and Whelan 197^) from spectra. 

We have explored this range of spectrcLL; types by employing the gravity darkening 

coefficients appropriate > to the endpoints 06.5 III-V and 09*5 III-V. The values 

! 

of (*^ 0 * are (0.295> - 0 .^ 5 )' aJid (0.1f60, ?-0.45) respectively. It is found 
that the gravity darkening coefficients eore essentially’ identical for luminosity 
classes III and V and .that the spectral: types near BO lb have gravity darkening 
coefficients within this reuige. The limb darkening coefficient has been taken 
to be 0.6 (Grygar 1965 ). | 

The best fit to all Ul observations for the above range of gravity . darkening 

I 

I / 

coefficients is shown as line segment a in figure 2. The small q endpoint is 

i ^ ^ 

the best fit 09*5 III-V solution and the large q. endpoint is the best fit O 6.5 

i . 

III-V solution.. The error ellipses on these endpoints are the formal errors 

'2 

from the x solution. The O 6.5 III-^V emd 09*5 III-V solutions represent the 
observations equally well, and . thus we cannot distinguish between the published 
spectreJ. types on the basis of the light,, curve analysis alone. 

Consideration of the residuals from -the beat -fitting synthetic light curve 
shows that the light curve provides a good representation of the observations 
since there are no systematic deviations. However, the internal errors do not 
provide a, good estimate of the , deviation of the observations from the best fit 
light curve. Thi^.,can be interpreted as evidence that the primary is irregularly 

varying in brightness with a dispersion of approximately 0.011 magnitude. 

i 

To test the sensitivity of the light curve parameters q and sin i to this 
variability* we fit the, 37. obse^rvationa obtained by deleting four discrepant 
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observations near the primeory and secondary minima, since the light curve 
parameters have a maximum sensitivity [to observations at the extrema* The 

' . i 

results of the fitting procedure for these obsLcrvatlons are shown as line 

' ^ 

segment b in figure 2. These results |lie outside the formal lo solution 
region based on all 1»1 observations, but not greatly so. 

I ‘ ; 

Thus, the present observations can be represented by the model developed / 
by KopaJL and Kitamura (1968) with a mass ratio near 0.123 and sin i near O.98U, 
the exact values depending upon the spectral type of the primary and the inter- 

: ' i ' 

pretation of the observations which are far from the mean light cturre. These 

two ambiguities lead to a 24. percent uncertainty in q. and 2 percent uncertainty 

i 

in sin 1. ' 

i 

,,i These values pf. q. and sin i are minimum values in that we have assumed 

^ ! 

\ 

that the Roche Ipbe is filled. Davidson and 0 striker (1973) have suggested 

I 4 

that early spepti;'?^ type, members of close binaries may rotate at less than the 

■ I 

I ' I 

synchronous velocity in which c^se it} is possible for a star to be in dynamical 
equilibrium, at sizes larger than the Roche lobe. The size mass ratio of 

a star principally affject the mean amplitude of the light variation, whereas 

i 

the relative depth of. the two minima is principally determined by the Inclin- 
ation .of the system. We therefore believe that our lack of knowledge about 
whether the Roche lobe is filled (or overfilled If the tidal lobe model is 
appropriate) p)[^puld not affect our estimation of the uncertainty in sin i. 


To estimate the uncertainty in q 


due to our lack. of knowledge of the 


size of Krzeminski*s star we note that the lowest order equations for the 
light changes of a close binfury (Kopai 1959* P* 412) can be used to find the 
following useful interpolation formula: . 

... « Bln^ 1 <1 (X *. T,j) ■ . 
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where the amplitude, Al, is given in terms of the radius of the star, r (in 
units of the separation); the limh darkening coefficient, x; and previously 
defined quantities. Since, near q *= 0.1 the tidal lobe radius is T percent 

. I 

greater than the Roche lobe radius (Davidson and Ostriker 1973) q may be 
smaller by as much as 21 percent. The approximate q limits for the tidal 
lobe model for solution tracks a and ^ are shown in Figure 2 as line segments 
^ and ^ respectively. Underfilling of the Roche lobe seems less likely since 
the spectroscopic observations of Osmer, Hiltnef, and Whelan (197^) constrain 
q from being much greater than 0.1. | 

We have the total imcertalnty of q to be 1*5 percent and of 

sin i to be 2 percent. It is expected that further observation and study of 
the system will j greatly reduce the present, uncertainty. 

I ' • ■ I - 

, I ( 

i IV. IDISCUSSION 

! 1 

; I 

■ I 

i 

The most important result which follows from the present study is an 
observational determination of the masses of Krzeminski's star and of Cen X-3 
and of. their separation. From the X-ray pulsation mass function and the 

projected size ol* t^e Cen X-3 orbit the following relations for the primary mass, 

; I 

M, and the ^separation of the coQiponents, a, are found: 

2) ' M = 1515 (1 + q)^ "/’ sin^ i 

j 

3) a = 1711 (1 + a) / sin i , 

i 

! . i 

where the mass and separation are in solar units. It can be seen that the 

i 

primary mass and the separation are principally sensitive to sin 1. 

. i. ... . . ij-. M . t j . , . , 

Squations 2) and 3) 5 the equivalent radius, R, of the Roche lobe (Kopal 1959) 

t I . !■ . t ( hi i . . I - . • : ; . , ■ • , 

/ 
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i 

in units of the component separation; and the radial velocity semi-amplitude t 

! • i 

K, of Krzemlnski's star In terms of the X-ray source velocity yield the results 
shown in table 1. 

j 

Table 1 


/ 


a « 19-5 i 1*2 Rj 
B » 12.8 t 0.8 R, 


K « 51 ± 19 km sec 


-1 


I q a 0.12 ± .05 i 

; sin i = 0.985 ± 0.02 j 

“primary " 20.5 ± 3.5 i 

j 

M = 2.5 ± 1.1 ! 

X 0 

The quoted uncertainty in each quantity reflects our estimate from the previous 

' (.t. i . I 

section of the Roche lobe filling uncertainty, the spectral type uncertainty, 

; • I 

and the observational interpretation uncertainty. 

• ,1 

It is interesting to compare these yG^.ues with stellar masses and radii 
tabulated as functions of spectral type. The masses of 06.5 III-V models are 
given as approximately 30 (Stothers 1972) and U5 (Robertson 1973). From 
equation 2) the value of sin i necessary to bring agreement between the model 
masses and the observed mass can be found. Lines and £ in figure 2 show 

i J 

these values of sin i for *^5 and 30 respectively. Clearly no agreement 
can exist with the present observations. However, the predicted masses of 

* t * * / » 1 1 . 

J 

of 09.5 III-V stars are in accord with the mass found here. Effects of mass 

t f ^ i /.,v ix-v-J.. -■ ' : .. ..... w 

1 

exchange on close binaiy members are not well understood and it is therefore 
unwise to rule out an atmospheric quantity, the spectral type, on the basis of 

a mass determination. We conclude that Kfzeminski's star cannot be an 06. 5 

• • ' • “ ‘ ' ■ i ■ • ■ ■ ‘ 

i 

star of "normal" mass. Radii of OB stars tabulated by Panagia (1973) show the 

, . I . . I • . I iii'i I • . 1 1 . ' ' 

radii of O9-BO supergiant a to be t%rice tl^t of KrzeminsM’a star. We, therefore. 


•Jin;. 



1 

i 

i 

i 

i 

i 

i 

I 

I 

i 

i 

! 

I 

I 

1 


j • - 9j- 

I I 

conclude that Krzeminski*s star cannot he a normal 09-B0 supergiant* Main 

’ . i 

sequence and giant stars in the range 06*5 - 09*5 are compatlhle vlth the 

! ^ 

' Hoche lohe diameter determined here. I 

i 

Prom the Doppler variation of the 4?8 X-ray pulsations it is found that 

1 ■ -.1 
the line of sight velocity variation of’Cen X-3 is Ul5 km sec (Schreier 

^ al * 1972). The line of sight velocity variation of the optical primary 

j ' , ' ' 

is simply q times this value. Though the observations of Osmer, Hiltner and 


Whelan (197^) do not place a strong constraint on the radial velocity variation 


of Krzeminski's steir, the radial velocity variation determined from the light 


/ 

/ 


curve analysis presented here is consistent with these, observations. 

^e duration of X-ray eclipse plaqes constraints on th^ size of the 

I 

‘ occ^ting object and the inclinjation of [the system. For purposes of compcnrison 

1 

ve may assume that the size of the occulting object is given by the Roche lobe. 

! ! 

The observed eclipse duration then constrains th^ possible values of q and sin i 


to a lociiB in the (q> sin i) , pJL^e. , For the duration of X-ray eclipse we have 
taken 0^U88 (Schreier et al. 197^), excluding the transition period. Cxirve ^ 


in figure 2 shows the :|*esult of a calculation in which, we have employed the 

J 

I 

third order expression for the surface radius given by Kopal (1959» P- 131). 

i 

I 

The mass .^'atio and incllni^t^pn found from the light curve do not fall on this 
locus. A possible resolution, pjT the discrepancy can be made by assuming that 
the X-rey op ^ulting surface is, larger th^ the Rwhe lohe hy approximately 13 
percent while maintaining the, photosphere, at the Roche lobe. The X-ray occulting 
surface would jnecessarily not be in d^amical equilibriiim as s\iggested previously 
for HD 77581 — 3U 0900-40 by Petro andj Hiltner (1974). An alternative explan- 
ation follows from the tidal lobe suggeLtion of Davidson and Ostriker (1973), 

i ■ ■ 

; • ^ . i I 1 > , . . i u . • • • • , 


I 
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namely that the X-ray occulting aiirface and the photosphere are closely coin- 

I 

cident €ind €u:e both near the tided lobe limit. 


V. COKCUJSION 

I 

I 

■ ■ I 

The lack of firm understanding of the structure of upper main sequence 

> I . I . . I .. .j I ... I . 

stars and of the evolution of these steucs in close binaries precludes the use 
of the results of this study to dlscrimlhate between the widely disparate 
spectral types suggested for Krzeminski’s star. The mass of . Krzeminski’s star 

is established! as 20,5 ± 3. 5 while that of Cen X-3 is estab-Lished as 2.5 ± 

' 1 '■ 

Mq. Resolutioni.of the,, spectral type ambiguity, further light cprve obser- 

i 

vations, and further radial velocity observations should lead to a much better 
determination of all parameters discussed in this paper. 
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FIGURE CAPTIONS 


Differential V light curve of Krzemlnski ' s star* Also 

i * 

shovn is the best fitting synthetic light curve for all 

. i 

4l observations* I 


Constraints on the system Cen X-3 - Krzeminski’s star, a 

1 

is the 06.5 III/V - 09.5|lH/V light curve solution track 

i 

for aJ.1 4l observations . j ^ is the same track for 3T of 

i '■ ■ 

the observations* The error ellipses shovn alb the endpoints 
depend only on the internetl errors of the obsearvations* 

t 

Lines £ €tnd d are . approximate lower limits to the value 
of q if the. tidal lobe model is correct. Locus : ogives 
the (q, sin i) values compatible with X-ray eclipse being 
caused by material at. the Roche lobe. Points to the right 


of this locus ! correspond 


to I eclipse being caused by . material 


beyond the, Roche lobe. Lines f and are the approximate 

j 

values of sin i necessary if the primary mass is to be 45 M. 


or 30 Mg respectively. 
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